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Abstract 

Mathematical models of quantum computers such as a multidimensional quan- 
^ I tum Turing machine and quantum circuits are described and its relations with 

lattice spin models are discussed. One of the main open problems one has to solve 
^^, if one wants to build a quantum computer is the decoherence due to the coupling 

with the environment. We propose a possible solution of this problem by using 
a control of parameters of the system. This proposal is based on the analysis of 
the spin-boson Hamiltonian performed in the stochastic limit approximation. 
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1 Introduction 

Quantum computation touches upon the foundations of mathematics, quantum physics 
and computer science. Various mathematical and physical aspects of quantum compu- 
tations have been discussed in many works, see for example [|3|-[0- Feynman [^ has 
pointed out that there is no efficient way of simulating a quantum mechanical system 
on computer and suggested that, perhaps, a computer based on quantum principles 
might be able to carry out the simulation efficiently. One hopes that remarkable prop- 
erties of quantum systems such as quantum parallelism and entanglement will lead to 
a major breakthrough in computer science. 

Quantum computers violate the modern form of the Church- Turing thesis. It was 
shown that certain problems, in particular the recursive Fourier sampling problem, 
can be solved in polynomial time on a quantum Turing machine, but relative to an 
oracle, requires superpolynomial time on a classical probabilistic Turing machine [Q. 
An efficient quantum algoriphm for factorizing integers was found by Shor ^j. 

There exist already small quantum computers working with few qubits. However it 
is not yet clear whether large quantum computers suitable for practical use are possible 
to build. One reason that quantum computer is difficult to build is decoherence of 
quantum superpositions due to the interaction with the environment. Decoherence 
can be defined as the decay of the off-diagonal matrix elements of the density operator 
in the computational basis. By using simple models like spin system coupling with 
quantum field it was found that decoherence is quite large. Different methods were 
discussed to reduce the decoherence in quantum computers ||TU[-|T^, still it is considered 
to be as one of the main open problems in this topic. 

In this work we propose an approach to the problem of reducing decoherence in 



quantum memory by using a control on the parameters of the system [15]. We show 



that one can choose the parameters of the system in such a way that decoherence is 



drastically reduced. This proposal is based on the previous analysis ||T6| of the spin- 
boson Hamiltonian performed in the stochastic limit approximation. Considerations of 
the decoherence problem in quantum computers performed in |T0|, ITl|, KM used a special 



case of the spin-boson interaction when no spin-flip transitions (or tunneling, in another 
interpretation) are present. Here we consider the complete spin-boson model, including 
the spin- flip term, and show, follow ||16[ , that this term plays in fact the crucial role for 
the reducing decoherence. Recently Viola and Lloyd |T^ have proposed to use the spin- 



flip transitions for the dynamical suppression of decoherence. The control procedure 
in their scheme is implemented as a sequence of radiofrequency pulses that repetitively 
flip the state of the system. 

Turing machines and circuits are two well known models of classical and quantum 
computers [Q. Models for the hardware of computers are given by circuits that compute 
Boolean functions. For a given sequence of Boolean functions /„ one looks for circuits 
Sn computing /„ with small size and small depth. A hardware problem is a Boolean 
function /„ and obviously each fn is computable. A software problem is a language 
L, i.e. a subset of the set of all flnite — 1-sequences and it is known that many 
languages are not computable. Circuits Sn for /„ are efficient (or uniform) if there is 



an efficient algoriphm for L, the union of all /^^(l). A circuit works only for inputs 
of a definite length whereas a reasonable program works for inputs of arbitrary length. 
Models for the software of computers are given by Turing machines. A universal TM 
models programmable computers. Uniform circuits can be simulated efficiently by 
Turing machines. 

In the next section we discuss mathematical models of quantum computers such as 
a multidimensional Turing machine and quantum circuits and its relations with spin 
and gauge theories on the lattice. Then the problem of decoherence is considered. 

2 Quantum Turing machine 

A quantum Turing machine (QTM) is obtained by quantization of a classical Turing 
machine (TM), so let us ffist describe the classical TM. If E is an alphabet, i.e. a finite 
set of simbols, then S* denotes the set of all strings (words) over E. Let S and Q be 
two different alphabets. 

A Turing machine M is a triple M = {Q, S, 6}. Here Q is a set of states, including 
the initial state go and the final state gi and S is the tape alphabet including a blank 
symbol Qq. The program is given by the transition function 

5 : (g - gi) X S ^ g X S X {/?, L, iV} 

The tape of TM is a sequence of cells, i.e. the one-dimensional lattice Z. Each cell of 
the tape contains a letter of S. At the beginning the input x = {xi, ..., Xn) is contained 
in the cells 0, ..., n — 1, all other cells contain the blank symbol ag- The central unit of 
the machine is at each point of time in one state g G g. The state at the beginning 
of computation if defined by go. The machine has a head which can read one cell and 
which can move from the cell i G Z to the cell i + 1 or i — 1. The head starts in the 
state go. If the machine is in state g and reads a and if S{q,a) = {q',a',(T), then the 
machine replaces the contents of the considered cell by a', the new state of machine is 
g', and the head moves one step in direction a (R= right, L= left, N= no move). The 
computation stops in the state gi. The result of computation can be read consecutively 
in the cells starting at until the ffist cell contains the letter qq. 

For input x let t{x) be the number of steps untill the machine stops and let s{x) 
be the number of different cells which are scanned by the head of the machine. The 
computational time t and the space s of the Turing machine are defined by 

t(n) = max{t{x)\\x\ = n} 

where |x| is the length of x and 

s{n) = 'max{s{x)\\x\ = n} 

A TM decides a language L if it computes 1 ii x E L and ii x ^ L. A configuration 
of TM is a vector 

K = {i,q,aj,,...aj^) 



where i E Z is the position of the head on the tape, q E Q is the current state and 
Qj e S is the contents of the cell j. Equivalently a configuration can be written as 

K = it,q,F) (1) 

where F : Z ^ S is a function with only a finite set of values different from Oq. A 
multidimensional tape TM has a d-dimensional tape Z"'. A rf-dimensional tape TM 
can move its head in 2d directions. 

Now let us perform quantization of the classical Turing machine. Let S and Q 
be two languages described above. The Hilbert space TC of states of quantum Tur- 
ing machine (QTM) is formed by complex vectors ip = {iPk},'^k G C indexed by 
configurations K of the classical TM with the scalar product 

K 

In 7i there is an orthonormal basis {e{K)}, {e{K),e{P)) = 6kp where e{K)K' = S§,. 
Because of (|l|) one can write e{K) = e(z, q, F). 

A quantum Turing machine M is a 4-tuple M = {Q,Ti,7i,U} where 7i is the 
Hilbert space described above and [/ is a unitary operator in 7i that satisfies the 
following condition. There exist three functions 

M^ : g X s X g X S ^ C, a = 0,±l 

such that one has the following locality condition 

(e(2, g, F), Ue{t', q', F')) = [6l+'u,{z) + 5\r'u^,{z) + 5>o(^)] \{ 4i? (2) 

for all configurations {i,q,F) and {i',q',F'). Here z = {q,F{i),q',F'{i)). The product 
Hj^j ^F(j) nonvanishes only in the case if functions F and F' are perhaps different at 
the point i. Here C is the field of complex numbers computable in the polynomial time. 
One can prove that it is enough to restrict ourself to the field of rational numbers. The 
fundamental role of the rational numbers is emphasized also in the p-adic approach 
to mathematical physics |18|. It would be very interesting to explore interrelations 



between these two approaches. 

Unitary evolution operator Ut is defined as the product of f/, Ut = U...U,t = 0, 1, .... 
Computations on QTM are performed as follows. Let i = and go corresponds to the 
input Fo and let e(0, go, -fo) be the initial vector. We denote ipt = f^te(0, go, Fq). If the 
vector ipt admits an expansion of the form 

^, = 5:c(F)e(0,gi,F) 

F 

where every vector e(0, gi, F) from the expansion includes the final state gi, and more- 
over the vector ipT- does not admit such an expansion for any t < t, then the computa- 
tion is considered the finished in time t and the vector ipt is the result of computation. 
In this case the probability to obtain the result F is |c(F)p/iV where N = Y, |c(F)p. 



The classical TM can be considered as a particular case of QTM when the functions 
Ua take values and 1. 

We have discussed the QTM whose tape is the one-dimensional lattice Z. The 
crucial property of the unitary evolution operator U is the locality condition (Q). A 
QTM provides perhaps not the best model for quantum computations. One could 
argue that models of quantum computations with a local Hamiltonian and not with 
the local evolution operator would be more realistic physically. 

To extend the efficiency of quantum computations we can look for a QTM with 
the multidimensional tape. It is known that a classical TM with the multidimensional 
tape is polynomially equivalent to a TM with the one- dimensional tape. From the 
other side properties of one- dimensional and multidimensional quantum spin systems 
are drastically different. Therefore it is worth to study a multidimensional QTM. 
We define a QTM M with the multidimensional tape by quantizing a TM with the 
multidimensional tape. 

Definition. A quantum Turing machine M with the multidimensional tape T/' is 
a 4-tuple M = {Q, S, 7i, U} where the Hilbert space 7i has the basis {e{i, q, F)} with 
i G Z"', q E Q and F : Z"' — > S. The locality condition reads 

(e(z, g, F), Ue{z', q', F')) = [ Y[ 6l^u.{z)] J] €!^ 

(T:\i—i'+a\<l j^i 

for all configurations (i,q,F) and {i',q',F'). Here z = {q,F{i),q',F'{i)). 

Quantum Turing machines with two- and three-dimensional tapes have, perhaps, a 
better chance for a practical realization of quantum computer than the one-dimensional 
models. 

Concerning the complexity of the Ising model it is known |T^ that finding a ground- 



state in the antiferromagnetic case is P - hard but in the ferromagnetic case the problem 
is NP - hard for general graphs and is in P for planar graphs. 

Evolution operators with certain locality conditions on the plane lattice are studied 
in the theory of quantum integrable models pO| . 



One has the following example of the QTM. Let Ti be the Hilbert space with the 
orthonormal basis e{i,a) where z G Z and a = 1,2, ...,n. Let U be an operator in 7i 
satisfying the conditions 

{e{i, a), Ue{i', a')) = 5i+i,i/A„„/ + Si^i^i^B^a' 

for all i,i',a,a'. Here A = {Aaa') and B = {Baa') are some n x n matrices. Then 
the operator U is unitary if and only if the matricies A and B satisfy the following 
conditions: 

AA* + BB* = /, A* A + B*B = I, AB* = 

A quantum computation is a sequence of unitary transformations. Any unitary 
matrix can be approximated by means of the product of unitary matrices of a simple 
form. Such a representation is the quantum analogue of the representation of a recursive 
function in terms of primitive functions. It is known that the set {e^'^*"^|7Z G Z}, where 



^ is a fixed irrational number, is dense on the unite circle. This can be interpreted as 
saying that the 1x1 matrix e^'^*^ is universal for the set of all unitary 1x1 matrices. 
A unitary d x c?-matrix [/ is of a simple form, if it has (after possible reodering) a 
block-diagonal form such that every block is a 2 x 2-matrix of rotations 

/ cos — sin 6 \ 
y sin 6* cos 6* j 

or it is a number e*^ with some 6. 

One denotes QMTq the subset of QTM with the evolution operator of the described 
form. 

There exists a classical algoriphm Q (i.e. a classical TM) that for a given d x d- 
matrix U and e > it computes in the time polynomial in d and log 1/e unitary 
matrices of the simple form Ui, ...,Un, where n depends polynomially on d, such that 

\\U-U,---Un\\<e 

Therefore to perform an arbitrary quantum computation one has to build a QTM 
performing the simple unitary transformations. There exists a universal QTM Ai 
that performs an approximate computation of simple unitary transformations with an 
arbitrary small error. Moreover the unitary evolution operator of QTM Ai includes 
only the block-matrices of rotations on the fixed angle. The Hilbert space C^ is called 
qubit. Any d-dimensional unitary matrix may be written as a product of 2d'^ — d unitary 
matrices {quantum gates), each of which acts only within a two-qubit subspace. 

Let be given two QTM Mi and M2 with the same tape alphabet S. One says that 
QTM M2 approximates QTM Mi if the following condition is satisfied. Let Mi for 
an input e(0, qo, F) computes the output ipi in time t. Then M2 for the same input 
computes the output ^2 in time polynomial in t and moreover 

11^1 -^2! I < e- 

A universal QTM should approximate any QTM. To this end one has to enumerate 
and encode all QTM in such a way that any QTM can be considered as an input for the 
universal QTM. Every QTM can be represented as a finite set {Q,T,,Ua{q,a,q',a')}. 
All such sets can be enumerated by using the tape alphabet {0, 1, gq}. 

There exists a QTM A1 (universal QTM) ^ that approximates any QTM . More- 
over there exists a universal QTM A4 belonging to the class QTMg, where the number 
6, < 6 < 11/2 is such that cos^ and sin 6' are rational numbers; in particular one can 
take cos6' = 3/5, sin^ = 4/5. 

Therefore for the practical realisation of quantum computers it is sufficient to im- 
plement rotations on two qubits subsystems of a quantum system. 

3 Quantum circuits 

Quantum circuits are quantum analogues of the classical circuits computing Boolean 
functions. A classical circuit can be represented as a directed acyclic graph. Similarly 



a quantum circuit is a sequence of unitary matrices of the special form associated with 
a (hyper)graph. 

Let us consider a quantum system with the Hilbert space 7i from / qubits, i.e. 
Ti = H'^\ where qubit H = C^. It is convenient to treat 7i as the set of functions 
on the set L = {1,2, ...,/} with values in H, i.e. 7i = H®^ = H^. If (7 is a subset of 
L, g G L, then one has H^ = H^(^H^^^ . Let ?7 be a unitary operator in if®'', r < / and 
the cardinality of the set (7 C L is r = |5f|. Let us denote f/^^^ the corresponding unitary 
operator in H^ and let Ag(f/) be its extension to a unitary operator in H^ obtained 
by tensoring to the identity, i.e. Kg{U) = U^^^ (g) id\{H^\3y Let U = {V, W, ...} be a 
finite set of unitary operators (basis of quantum gates) each of which acts in one of the 
spaces H^", Va < I, a = 1, ...,q. 

A pair (L, F) where L is a set and F is a family of its subsets is called a hypergraph. 
If one of the members of the family is marked then such a hypergraph is called a 
network or circuit. We will deal with the hypergraph [L,r), where L = {1, ...,/} and 
F = {Li, ..., Lt} is a family of subsets from L, Li C L. 

Definition. A quantum curcuit Si is a 5-tuple 

Si = {{L,T), if^ U, Ur} 

Here {L,T), H^ and U are the described above hypergraph, Hilbert space and the 
basis of gates, and Ur is a unitary operator in if ^ of the form 

f/r = AL^(f/T)...AL,(f/i), 

where f/, are unitary operators from the basis U. It is assumed that the hypergraph and 
the basis are consistent in the sence that Ui acts in H'^°'i^) and \Li\ = ra(i), ^ = 1, ...,T. 
Here a is a function from {1, ..., T} to {1, ..., g}. 

If {ca}, a = 0, 1 is a basis in qubit C^ then the unitary operation in C^ of the 
form Ca ® Cb — >■ Ca ® ea0b is called the quantum controlled NOT gate. Here © denotes 
addition modulo 2. One can show that the quantum controlled NOT gate, together 
with one-qubit gates, is sufficient for any arbitrary quantum computation. 

By using the Dirac notations one writes |ai,...,a; >= Cq^ ... Cq;. The basis 
{\a >= \ai,...,ai >} in if' is called the computational basis. Here a = I]'lo2*aj is the 
binary decomposition on the number a. 

More details about mathematical models of quantum computers one can find for 



example in [21 



4 Reducing of decoherence in quantum computers 

The maintenance of quantum coherence is a crucial requirement of the ability of quan- 
tum computers to be more efficient in certain problems than classical computers. One 
can simulate the environment as a classical or quantum white noise [^. In a simple 



model of spin coupling with the massless quantum field it was found |T0[ that for quan 



tum computations not only the coupling with the environment must be small, but the 



time taken in the quantum calculation must be less than the thermal time scale 1/T 
where T is the temperature of the initial state of the environment. 

The tape (memory) cells are taken to be two-level systems (qubits), with each of 
the levels having the same energy and the two states are taken to be the eigenstates of 
the spin operator a^- 

Any (pure or mixed) state of a quantum circuit Si can be described by a density 
operator of the form 

Ps{t)= E Pa6(t)|a><6| (3) 

a,fe=0 

where {|a >} is the computational basis in if'. The degree of the quantum coherence 
is described by the off-diagonal elements pah-, a- 7^ h oi the density operator. The 
decoherence is the decay of the off-diagonal elements of the density operator in the 
computational basis. 

In the simple model of the quantum computer interacting with the environment 
(reservoir) represented by a quantum field (a family of harmonic oscillators) one as- 
sumes that the total Hilbert space y^ H^ ® T where H^ is the /-qubit space while T 
is the bosonic Fock space. If pTot{t) is the density operator of the total system then 
to get the density operator ps(t) of the quantum computer one has to take the partial 
trace over the reservoir space 

Ps{t)=TrApTot{t)) 

The Hamiltonian decsribing the coupling of qubit with the environment (the spin- 
boson Hamiltonian) has the form 

Hx = -^Aa, + l^^z + J dku{k)a+{k)a{k) + XaMiol + ^^(^)) (4) 

where ax and a^ are Pauli matrices, e and A are real parameters interpreted respectively 
as the energy of the spin and the spin-flip parameter. Here 

A+ig) = j a+{k)g{k)dk, A{g*) = J a{k)g*{k)dk 

where a{k), and a~^{k) are bosonic annihilation and creation operators 

[a{k),a+{k')] =6{k-k') 

which describe the environment. 

The one-particle energy of the environment is denoted uj{k) and we assume u{k) > 
0. The function g{k) is a form factor describing the interaction of the system with the 
environment, A is the coupling constant. It is well known that, in times of order t/A^, 
the interaction produces effects of order t. Thus A provides a natural time scale for the 
observable effects of the interaction system-environment. 

Leggett et al. [Q have found a very rich behavior of the dynamics of the Hamilto- 



nian (^) ranging from undamped oscillations, to exponential relaxation, to power-law 
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types of behavior and to total localization. Main qualitative features of the system 
dynamics can be described in terms of the temperature (i.e. the initial state of the 
environment) and of the behavior, for low frequencies cu, of the spectral function 

J{u) = f dk\g{k)\^6{uj{k) - u) (5) 

The dynamics of the Hamiltonian (^ in the so called stochastic approximation have 
been considered in [0.1t was found that the pure oscillating regime when no damping 
and no decoherence is present is descibed by the simple equation 

J(z/A) = (6) 

where 

ex2 



.= W1+,^ 



The basic idea of the stochastic approximation (see [Q) is the following. If one 
has a Hamiltonian of the form 

Hx = Ho + XV (7) 

then the stochastic limit of the evolution operator 



U^^\t) 



^UHo^-itH, 



A 



is the following limit (when it exists in the sense of the convergence of matrix elements): 

The limiting evolution operator U{t) describes the behavior of the model in the time 
scale t/X^. 

In order to apply the stochastic approximation to the Hamiltonian (|^), we write 
(|^) in the form (^ where 

Hq = Hs + Hr 



The system Hamiltonian Hg is 



Hs = -ijl^^x + 2 ^^2 



and the reservoir Hamiltonian Hb. is 

Er= f dkuj{k)a+{k)a{k) 
The evolution operator U^^\t) satisfies equation 

'-^ = -^xv(t)u^^Ht) 

at 



where V{t) = e^^^°Ve ^*^° has the form 

V{t) = a,{t){Aie-''-g*) + A+{e''-g)) 

and 

a,{t) = e''^^a,e-''"^ (8) 

Let us compute (^ . The eigenvahies of the Hamihonian Hs are 

Hs\e± >= \±\e± > 



where 



and 



1 , 1/1 

A± = ±- Ai^, e± >= : 

2 1^ Jl + J4^^'T 



e , /. /e^2 



Notice that: 



^±=A^"' " V* ' VA 



1 -/i^ 

e±\(T^\e±) = ^ ^ J ; (e+|(T^|e_) = (e_|(T2|e+) = 1/u 



l+Z^T 
Therefore 

^Jt) = LU^ DD+ + ^—^ D+D + u-'e''"'' D + u-^e-'"'''D+ 
^' i + nl. l + /ii 

where 

D = |e+ >< e_| 

The interaction Hamihonian can now be written in the form : 

3 

V{t) = ^(Z}+ ® A(e-'*""(7*) + h.c.) 

a=l 

where the three spectral frequencies correspond respectively to the down, zero, and up 
transitions of the 2-level system, i.e. 

uJi{k) = uj{k) — v^ , uj2{k) = uj{k) , uj^{k) = uj{k) + v^ 

D, = u-'D+ , D, = \^ DD+ + \^D^D , D, = v-'D^ 

1 + /il 1 + Ai+ 

The spectral frequencies uJ2{k) and uj^{k) are positive. Therefore in the stochastic limit 
one gets only one white noise field. Still a remnant of the interaction remains because, 
after the limit, the system evolves with a new hamiltonian, equal to the old one plus a 
shift term depending on the interaction and on the initial state of the field. This was 
called in [jl6| a Cheshire Cat effect. 
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The limiting evolution equation can then be written: 
dU{t) 



dt 
where 



Dh+{t)U{t) - D+U{t)b{t) - (7 + ia)D+DU{t) - iipU{t) (9) 

7 = Z/^^TT J(z/A) , 



,2 N 2N 



.^.-(/(-.A)-/(.A)).||i^) -(i^l 1/(0). 

and we denote 

diu'J{iu') 



dk\g{k)fSiujik)-uj) ; I{uj) = P.P. 



Cij' — CiJ 



where P.P. means the principal part of the integral. The operators 6(t), h^{t) satisfy 
the quantum white noise relations 

In the notations of quantum stochastic equations the equation (|^) reads 
dU{t) = {DdBt - D+dBt - (7 + ia)D+D - i^)U{t) 



Notice that all parameters 7, a and ip in the evolution equation ( pTD are expressed in 
terms of the spectral density J{oj) and parameters A and e of the original Hamiltonian. 
For zero temperature the stochastic approximation to the vacuum expectation value 
of the Heisenberg evolution of o"^ is given by 

P{t) = {U\t)aMU{t)) 

From equation (^ one gets the Langevin equation for P{t) which solution is 

P{t) = i^-^e-^'^D+e'^"-^^^' + Z}e-^("~''^)*) + 

We obtain the pure oscillating behaviour and no decoherence if 

7 = z/~VJ(z/A) = (11) 

For a non-zero temperature we get a stochastic evolution equation of the same form 
as before only with new constants 7, a and ip. More precisely: 

2 .._ .^__.. /?i^A 



7 = z/ 7r J(z/A) coth 
11 



a 



'i-^^lV A--"^ 



(/+(0) + /_(0)) + 



+z/-2(/+(-z/A) - /+(z/A) + /_(-z/A) - /_(z/A)) 
where spectral densities are 

Here J{co) is the spectral density (j^) and /5 is the inverse temperature. The functions 
I±(^co) are defined by 

Uu) = P.P. f '"^-^^("^^ 
J co' — CO 

One has the same as for the zero-temperature expression (0) for P{t) but now with 
new constants 7 and a depending on temperature. The condition for the reducing of 
coherence is still the same (0). It seems this condition is a rather week requirement 
on the parameters of the interaction between quantum computer and the environment, 
so one can hope to use it to reduce decoherence. 
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